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Engineering applications for aircraft noise prediction contain models for physical phe- 
nomenon that enable solutions to be computed quickly. These models contain parameters 
that have an uncertainty not accounted for in the solution. To include uncertainty in the 
solution, nondeterministic computational methods are applied. Using prediction models for 
supersonic jet broadband shock-associated noise, fixed model parameters are replaced by 
probability distributions to illustrate one of these methods. The results show the impact of 
using nondeterministic parameters both on estimating the model output uncertainty and 
on the model spectral level prediction. In addition, a global sensitivity analysis is used 
to determine the influence of the model parameters on the output, and to identify the 
parameters with the least influence on model output. 


Nomenclature 


A Area 

B Intercept parameter, see equation (3) 

D Diameter 

E Amplitude power exponent parameter, see equation (6) 
Jo Bessel function of the first kind, order zero 

Ji Bessel function of the first kind, order one 

L Gaussian source half-width function, see equation (7) 

M Mach number 

N Total number of random samples 

P Pressure 

R Radial distance 

Sti Total effects sensitivity measure, see equation (1) 

St Strouhal number 

T Temperature 

U Velocity 

Xi 1-th parameter in vector X 

Y output quantity from model with parameters X 

X Vector of parameters 

X* X fixed at a set of values 

E{ } Expected value of the enclosed quantity 

Var( ) Variance of the enclosed quantity 
2 

A Shock cell strength amplitude function, see equation (7) 
c Amplitude factor, see equation (7) 

c Speed of sound 
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/ Frequency 

Jo,;/ /-'-th zero of a zeroth order Bessel function 

to Shock power factor 

(p 2 ) Mean-square acoustic pressure 

x Independent regression variable 

y Dependent regression variable 

z Linear regression residual error, see equation (5) 

Subscripts 

oo Ambient conditions 

i/, N v integer index, maximum value N„ 

All parameters in vector not including the *-th parameter 
c Convective velocity or Mach number 

d Design conditions 

e Nozzle exit conditions 

i,I integer index, maximum value I 
j Fully expanded jet conditions 

k , K integer index, maximum value K 
r Randomly selected variable 

sh Shock-associated noise related 

t Total conditions 

Symbols 

(3 Shock parameter 

Xn Chi-squared distribution with n-degrees of freedom 
e Nonlinear regression residual error, see equation (6) 

7 Ratio of specific heats 

pi Density with one subscript 

Pi 2 Correlation coefficient with two subscripts 

a Standard deviation 

6 Inlet directivity angle 

N(p,a 2 ) Normal distribution with mean p and variance a 2 

I. Introduction 

The development of tools to computationally predict the sound generated by fluid flow or when fluid 
interacts with a surface continues to be an important part of the research effort for developing future envi- 
ronmentally friendly aircraft. For aircraft propulsion system related jet noise, a wide range of computer-based 
schemes are included in a hierarchy of these noise prediction tools from numerically intensive, extremely long 
running, direct simulations based on solving first principle equations to near instantly computed, empirical 
methods using experimental databases. The numerically intensive computations may simulate all the proper 
physics for jet noise predictions, but running them is so costly that they currently have little practical usage. 
At the other extreme, the empirical models are limited in their ability to predict jet noise beyond a relatively 
small set of geometries and operating conditions with which they were developed. In between these two ex- 
tremes are found methods that use some combination of numerical techniques and physics-based models in 
an attempt to increase the practical ability to predict jet noise over a wider range of jet operating conditions 
and nozzle geometries. Among these are included tools to compute the unsteady fluid flow using various 
types of filtering and turbulence modeling schemes to obtain a noise source field that is used to compute the 
radiated noise field, and tools that use acoustic analogy source models based on mean fluid turbulent statis- 
tical values to compute a solution for the radiated noise field. Even though these latter numerical methods 
still include some level of modeling that reduces total computation time, the time required to compute the 
unsteady or the steady mean flow can inhibit the number of cases run for a design study. 

For engineering applications in aircraft noise prediction at the component and especially at the system 
level, extensive modeling of the physical phenomena is required to obtain many solutions for design, de- 
velopment, and parametric studies. Analytical development and intuition are used to develop models that 
describe physical phenomena. In the end, these models will contain parameters that must be determined 
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by optimization and calibration. These parameters are found by using numerical simulation and/or exper- 
imental data. Regression is often used to set these parameters. This deterministic approach may provide 
reasonable solutions, but usually within limited conditions. However, the determination or the calibration 
of these parameters is often inexact especially when using experimental data that contains uncertainty, that 
may be insufficient for good statistics, or that may be of poor quality. Thus there is a question if the best-fit 
parameters are truly the optimal parameters for the model. It may be more appropriate to include the pa- 
rameter uncertainty in the model 1 and/or to compute the solution using nondeterministic methods. 2 In the 
present paper, we consider the introduction of nondeterministic parameters into models for the prediction 
of broadband, shock-associated noise from supersonic jets. This was done as part of a broader study and 
assessment of jet noise prediction models both as a method of determining the quality of the prediction and 
as a method of providing better prediction capability. Nondeterministic approaches are shown to impact the 
predicted results and can be used to estimate the output uncertainty and to determine the influence of the 
model parameters on the output. 

A brief background for this paper is given in the next section. A previous assessment study of jet noise 
prediction models, including models for broadband shock-associated noise, determined the uncertainty in the 
output from the prediction models. The shock-associated noise prediction models were easily modified to 
apply the method of computing with nondeterministic parameters. This is followed by a section describing 
some of the details of the prediction models. One model relies on a two-parameter regression analysis of 
experimental data to provide a scheme to predict noise amplitude as a function of jet operating conditions 
and directivity angle. The regression analysis provides the statistical quantities to define the distribution for 
the uncertain parameters. The other model has a large number of parameters with unknown distributions. 
The results are then presented in three parts. The first discusses the effects of the input operating condition 
uncertainty on model output uncertainty. Next, the results for the model using the regression analysis is 
presented. The main issue is which type of regression is appropriate. The last part presents results for the 
many parameter model. This part focuses on the choice of distribution for the uncertain parameters and the 
determination of the influence of the parameters on the output. This helps identify parameters that most 
affect the output uncertainty. 


II. Background and Methodology 

Bridges et al. 3 assessed the capability to predict jet noise and showed results from empirical and acoustic 
analogy based tools compared to measured data over a wide range of operating conditions. A general 
conclusion was that empirical tools could make overall sound pressure level predictions to within ±6 dB 
of measured values for a broad range of nozzle shapes and operating conditions when averaged over all 
directivity angles. For simple axisymmetric nozzles, the predictions improved to agree within ±1.2 dB. In 
comparison, an acoustic analogy tool was only applicable to a limited range of conditions consisting of 
axisymmetric, cold subsonic jets. The accuracy for these comparisons were also within ±1.2 dB. That study 
showed some of the issues related to how the results from different jet noise prediction models compare 
to the same set of observations and how to judge the quality of those predictions. Khavaran & Bridges 4 
further analyzed the measured jet noise data at supersonic conditions, separating mixing noise and shock- 
associated noise spectral components from the total jet noise spectrum, to provide better data for jet noise 
model predictions and improvements. The results of that component separation allow us to focus this 
study on shock-associated noise prediction capabilities for jets from convergent nozzles. To illustrate the 
results from different models, example spectra that compare predictions for shock-associated noise with the 
shock-associated noise component of measured data are shown in Figure 1. The plots show power spectral 
density as a function of the Strouhal number for under-expanded jets with a fully expanded Mach number 
of about 1.5. There are two jet temperatures represented at two directivity angles from the inlet of 50 and 
90 degrees. Three predictions are shown from the following models: The analytically based shock-associated 
noise prediction model from Tam, 5 the intensity scaling prediction model based on the work of Khavaran & 
Bridges 6 and further detailed in Appendix A of this paper (labeled sJet ), and the Reynolds-average-Navier- 
Stokes-based (RANS) model of Morris & Miller. ' Three different models attempting to predict the same 
experimental results. 

The shock-associated noise models under consideration here were originally developed to make determin- 
istic, one-time calculations using parameters found through least-squares fitting, averaging, point matching, 
or other methods to calibrate computed predictions using experimental data. These fixed parameters in 
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reality have an uncertainty associated with them that is ignored in the basic usage of these prediction mod- 
els. To incorporate the uncertainty, knowledge about the quantities that describe the uncertainty must be 
available and/or the distribution of the uncertain parameter must be known or estimated. This informa- 
tion may be unknown or hard to obtain. In the case of the sJet code, the uncertain parameters are well 
described in terms of their probability distributions and the moments associated with those distributions. 
Thus, it is readily possible to replace fixed parameters with their nondeterministic distributions. Once this 
is done, a method must be applied to perform the probabilistic analysis. One approach is to perform Monte 
Carlo sampling using the original prediction model. Depending on the scale of the problem, the number 
of computations may become quite large. Latin hypercube sampling helps reduce the number of sampling 
computations required to achieve convergence for good statistics. 8 Alternatives are to replace the original 
computational model with an analytical model that is more computationally efficient. Such methods include 
first and second order reliability methods, advanced mean value methods, and the response surface method. 9 
The analytical model may be faster to compute, but it is approximate and contains inaccuracies. The two 
models to be used here, the sJet and Tam models, are already approximate and are easily amenable to 
replacing fixed parameters with nondeterministic distributions. The RANS-based model was not available 
for this study. Thus, the Monte Carlo approach was chosen. A software tool that easily performs the Monte 
Carlo calculations using the external jet noise prediction codes is called DAKOTA. 8 This program can be 
set up to perform the following: 

1. Generate random inputs from predefined distributions. 


2. Execute the prediction model code using the random inputs. 

3. Perform statistical analysis on the outputs from the model code. 


In addition to this Monte Carlo calculation, the DAKOTA code can perform sensitivity analyses of the 
parameters to determine their influence on the output, optimization studies to find the best parameters, and 
other related computations. 

Sensitivity analysis is often used in conjunction with uncertainty analysis. While uncertainty analysis 
attempts to quantify the level of model output uncertainty given a set of model parameter uncertainties, 
sensitivity analysis is concerned with attempting to determine how much influence a particular parameter 
has on the output or output uncertainty and then rank the parameters in degrees of importance. This allows 
parameters of lower importance to be fixed and more attention and resources can be focused on parameters of 
higher importance. One sensitivity analysis method supported by the DAKOTA code is called variance-based 
decomposition. Let the parameters be denoted by the vector A". If a particular parameter Xj was absolutely 
non-influential, then the variance of the output Y due to changes in X, when all the other parameters in 
X are fixed, denoted as XX, must be zero. We compute the conditional variance Varx^Yj-X^i = XX) 
where the subscript X, indicates the one parameter in X that is being randomly sampled while all the other 
parameters are held fixed. The conditional variance must be zero for a second set of XX, a third set of 
XX, and so forth. Averaging these conditional variances over the set X^i will not change this result, hence 
the average will be zero, also. Comparing the average to the total variance of the output Y , we obtain the 
total effect sensitivity measure Sti for the parameter Xi, 


Sri = 


Ex^,{Var, Y ,(F|A^)} 

Var(K) 


(1) 


where Sti includes both a measure of the direct effect of Xi and a measure of the interaction effects of 
Xi with all the other X^i on the output variance. In actual calculation, Sxi will most likely be nonzero; 
hence, a small value for Sti represents an Xi parameter with less influence on the output than a different 
Xi parameter with a larger Sxi value. Since the computation of the total effect sensitivity measure involves 
varying all the parameters simultaneously, it is considered to be a global measure in contrast to a local 
measure that involves varying only one parameter at a time. The details of variance-based decomposition 
and its implementation to obtain sensitivity measures are found in Saltelli et al. 10,11 and the references 
therein. 


III. Models and Data 

This section introduces the two jet broadband shock-associated noise prediction models. The development 
approach for these two models were entirely different. The sJet code was developed from a study of the data, 
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particularly far- field spectral densities at various directivity angles. The Tam model was developed from 
derivations of the equations of motion. Only after this development were simplifications and adjustments 
made when results were compared to measured data. The section concludes with the experimental conditions 
for the data used for comparison with predictions in this paper. 

III. A. sJet Model 

In the jet noise prediction modeling approach proposed by Khavaran & Bridges, 6 scaling laws are used 
to derive relationships between fundamental jet parameters and the total radiated power computed by 
integrating over the appropriate power spectral density. To enable this process for imperfectly expanded 
supersonic jets, the radiated broadband noise was separated into distinct spectral components associated 
with mixing noise and shock-associated noise. Here, we are only concerned with the model for shock- 
associated noise. The total shock-noise intensity scales with the shock parameter /3, a measure of the shock 
strength in the flow given by the amount to which the fully expanded Mach number deviates from the nozzle 
design Mach number. In general, this may be expressed as 

( P 2 )sh OC (3 m where f3 2 = | M 2 - Mj | . (2) 

Based on experimental measurements, the shock power factor m varies with the noise radiation directivity 
angle and the jet temperature. 

Following Khavaran & Bridges, 6 a linear regression is applied in the form 


y i =m(e,T)x i + B{e,T) (3) 

where the dependent variable y,; and the independent variable Xi are given by 

Vi = OASPL s /,(0,T) and x t = lOlog/3,, *=1,2, ...,/. (4) 

The overall sound pressure level (OASPL) is for the shock-associated noise computed from the shock- 
associated noise spectral component of the total radiated noise over the Strouhal number range 0.1 < St < 10 
where St = fD e /Uj. The intercept parameter B is a measure of the shock noise overall level in decibels 
when the shock parameter (3 equals one. The results of the linear regression are the mean values rfi and B, 
the standard deviations a m and erg, and the correlation coefficient p m B between the m and B parameters. 
Once m and B are determined from multiple linear regressions covering a range of inlet directivity angles 6 
and total temperatures T t , the shock-associated noise levels are now set in general and it remains to compute 
the spectral shape. See Appendix A for the detailed equations. 

The linear regression model given in equation (3) is based on the general model regression equation 


Ui = m,Xi + B + Zi (5) 

where the Z{ s are assumed to be independent, normally-distributed, random values with zero mean and 
equal variances representing the random or residual error in the linear regression. 12 This assumption is most 
likely not strictly maintained in the above linear regression. Equation (3) is a logarithmic transformation of 
equation (2) . While this creates a linear equation from a nonlinear equation, the errors are also transformed 
and may violate the assumption for Z{ in equation (5). 13 In the original formulation, (p 2 ) s h is an overall power 
level computed by adding together the power levels in each frequency bin of the shock-associated noise power 
spectral density (PSD). Each power level at each frequency in the PSD has a chi-squared distribution y 2 
with n-degrees of freedom after spectral averaging. 14 A summation over K frequencies in a power spectral 
density with each power level having a y 2 distribution results in a total power level with a chi-squared 
distribution Xkxu with (K x n)-degrees of freedom. A nonlinear regression equation for the total power 
level can then be written as 

y i = 10^(0, T ) f 3 T ( 0 ,T ) + ei ( 6 ) 

where both y, = ( p 2 ) s h(0,T ) and e*, the residual error, have in general a yjf X ^distribution. Using 10-® is a 
recognition that (p 2 ) s h has a potential range covering many orders of magnitude. It can be shown that when 
weighting is used and the chi-squared distribution has a large number of degrees of freedom, the nonlinear 
regression formulation of equation (6) is appropriate for the given data and that the residual errors more 
closely have the desired normally-distributed properties. See Appendix B for details. 
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III.B. Tam Model 


The broadband shock-associated noise model of Tam 5, 15 is based on the theory that large-scale turbulent 
structures generated in the supersonic jet mixing layer interact with the periodic-like shock cell structure 
that exists in imperfectly expanded jets. This interaction causes pressure perturbations that radiate as 
sound. The shock cell structure is modeled as a wave guide containing a series of modes and the large 
structures are modeled as a sum of instability waves each with their own amplitude envelope and phase 
speed characteristics. These characteristics affect the amplitude and directivity of the radiated noise. The 
model is applicable to both under-expanded and over-expanded jets, but in this paper, it is applied only to 
under-expanded jets issuing from convergent nozzles. 

Following extensive algebraic manipulations and modeling approximations, Tam 5 derived an equation for 
the far field mean-square acoustic pressure 


(P 2 )k = 


Pj 


1+ 


1 -1 ^ 


— r2 , 


c(c b i) L 2 (h sl ,x c i,x c2 ) Aj A {a b i,a b2 )p 2 ao cl 0 U J 


2 


R 2 St k fk 


JiOV) 


exp [.Tv] (7) 


and 

T v 


fvjVcl, Vc2, Wn, W 12 , W 2 1, Wja) 
.fk 


- 1 


L 2 (h s i,x c i,x c2 ) 
2 In 2 


Uj 


Uc{Vcl,V c2 ) 


2 

(1 + M c (v c i,v c2 )cos9)' 2 


where pj, Uj, Mj , and Aj are the density, velocity, Mach number, and cross-sectional area, respectively, of 
the jet at the fully expanded condition, p^ is the ambient density, c 0 0 is the ambient speed of sound, R 
is the radial observer distance, 6 is the inlet directivity angle, fk is the fc-th frequency in a spectrum, Stk 
is the Strouhal number fkDj/Uj, Dj is the diameter corresponding to Aj, and f v is the peak frequency 
corresponding to the mode number v.. The denominator factor Stk in equation (7) is fixed at 0.35 when 
Stk < 0.35. In the summation of waveguide modes, Ji is a first order Bessel function and jo,v is the z/-tli 
zero of the zeroth order Bessel function Jo- The functional arguments with two subscripts represent the 
twelve parameters in equation (7) that were derived by fitting equations to experimental or calculated data. 
Table 1 defines these parameters and the detailed equations are given in Appendix C. 


Parameter 

Contributes to 

Model or Mean Value 

a b l 

amplitude 

3 (Mj/M d > 1), 6 (Mj/M d < 1) 

a b2 

amplitude 

3 (Mj/M d > 1), 5 {Mj/Md < 1) 

Wn 

1st waveguide mode wave number 

1.596 

w 12 

1st waveguide mode wave number 

0.1773 

W 2 1 

2nd waveguide mode wave number 

1.265 

W u 1 

v-th. waveguide mode wave number 

1.1 

Cbl 

amplitude 

2.65xl0" 4 

h s i 

source half-width length scale 

3.3 

X C 1 

potential core length 

4.3 

X C 2 

potential core length 

1.2 

V C 1 

convective velocity 

0.7 

V C 2 

convective velocity 

0.025 


Table 1: Tam model parameter values determined from experimental data or numerical calculations. 


III.C. Experimental Data 

The results of exploring the nondeterministic parameters in the shock-associated noise prediction models are 
compared to measurements consisting of noise data from supersonic jets issuing from a 2-inch exit diameter, 
circular, convergent nozzle. This nozzle was used in an extensive test program to measure the noise from jets 
at a variety of nozzle pressure and temperature ratios. 4 The measured noise spectrum from this convergent 
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Case 

Tt/Too 

Pt/Poo 

Mj 

Uj / Coo 

P 

7 j 

Too(K) 

-Poo (Pa) 

Coo (m/s) 

Cold 

1.01 

3.636 

1.493 

1.245 

1.109 

1.4000 

279.2 

97790 

334.93 

Hot 

2.16 

3.629 

1.494 

1.828 

1.109 

1.3942 

279.6 

97839 

335.18 


Table 2: Test conditions for a convergent nozzle used in the prediction models. 


nozzle at a variety of operating conditions were recently compared to predictions from both jet mixing noise 
and jet shock-associated noise models. 16 A subset of those conditions, one a cold jet and one a hot jet, are 
used in this study and are given in Table 2. Example spectral comparisons between measured and predicted 
shock-associated noise results at the operating conditions of Table 2 are shown in Figure 1. 

IV. Results 

The results of this paper are divided into three parts. The first part discusses the effects of jet operating 
condition uncertainty on the predicted results. This is followed by a presentation of results where the 
uncertain parameters have known distributions. In the final part, the prediction model has many parameters 
with no clearly defined distributions. Estimates are made for these distributions in order to complete the 
analysis. 

IV. A. Operating Condition Input Parameter Uncertainty 

The validation of a prediction code is performed by comparing computed values with the same values obtained 
from an experimental measurement. For jet noise prediction, this comparison is done for various noise 
level metrics including overall and spectral levels at specific jet operating conditions. Hence, the operating 
conditions from the experiments are used as input to the jet noise prediction code. The primary experimental 
measurements are pressures, temperatures, and lengths or distances associated with the geometry of the 
nozzle and the location of measurement microphones. Any code specific operating condition input parameter, 
such as Mach number, can be computed using the measured values. The purpose here is to demonstrate 
the uncertainty of predicted noise levels given the known uncertainty in the experimental jet operating 
conditions. For this example, we will use the Tam model. 

The ability to accurately measure pressures, temperatures, and lengths at the NASA Glenn Research 
Center aeroacoustic facilities has been documented in a recent assessment of noise prediction capability. 17 
Given those results, Table 3 shows an example of the measured operating conditions and the standard 
deviation a for the assumed normally distributed measured quantities. The table shows the measured total 
pressure and temperature, P± and T t , the measured ambient pressure and temperature, Poo and Too, the 
measured nozzle exit diameter D e , and the radial distance R from the center of the nozzle exit area to the 
microphone. The example conditions are for the noise radiated from a jet issuing from a convergent nozzle 
with the hot jet operating conditions listed in Table 2. 


Parameter 

Measured Value (mean) 

<7 

Pt 

355058 Pa 

139.3 

T t 

605.3 K 

1.89 

Poo 

97839 Pa 

139.3 

Too 

279.6 K 

1.89 

D e 

0.0508 m 

0.000051 

R 

5.08 m 

0.0127 


Table 3: Example measured jet operating conditions with standard deviations. 

Following the Monte Carlo procedure given above, sampling calculations were performed from 100 to 
10000 samples. In each case, the means and 95% confidence intervals were computed for the predicted 
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shock-noise spectra at directivity angles from 50 to 90 degrees from the inlet direction. From these results, 
the maximum width and the average width confidence intervals were determined across all computed spectral 
frequencies and directivity angles and presented in Table 4. The increasing number of samples N results 
in both more clearly defining the distribution of the computed solution and reducing the estimate of the 
standard error a/y/N, thereby reducing the confidence interval. The table clearly indicates that the effects 
of input parameter uncertainty on the predicted shock-associated noise, sound-pressure-level spectrum for 
the jet operating conditions listed in Table 3 are less than 0.2 dB at 95% confidence at the maximum and on 
average, less than 0.04 dB. These levels of uncertainty would not be noticeable on spectral density plots such 
as those shown in Figure 1. In general, this shows that the ability to make accurate measurements of input 
parameters will lead to small uncertainties in the predicted values. This does not mean the predictions are 
accurate as can be seen in Figure 1 where many decibels of error can separate the predicted values from the 
measured values. Thus, we cannot rely on determining prediction uncertainty solely from the uncertainty in 
the jet operating condition input parameters, but must consider the uncertainty of the parameters within 
the model. We turn to these parameters, next. 


No. of Samples N 

Max. 95% Cl Width 

Avg. 95% Cl Width 

100 

0.198 dB 

0.036 dB 

200 

0.146 dB 

0.025 dB 

400 

0.099 dB 

0.018 dB 

1000 

0.063 dB 

0.011 dB 

10000 

0.020 dB 

0.004 dB 


Table 4: Maximum and average widths of spectral 95% confidence intervals (Cl) at upstream angles. 


IV. B. Uncertain Parameters with Known Distributions 

The sJet code computes the shock-associated noise spectrum by first computing an overall level value and 
then combining that value with a spectral function to produce the final results. For the convergent nozzle 
used here, the jet operating conditions and directivity angle set the spectral shape. The only uncertain 
values are m and B that were determined by regression and only enter the model prediction through the the 
overall level in the first two terms on the right hand side of equation (16). Hence, we compare the outputs 
using the uncertain parameters m and B in those two terms to the output obtained using the mean values 
m and B. The spectrum result that would follow from equation (22) will have exactly the same level of 
uncertainty at all frequencies. 

The linear regression analysis performed in Khavaran & Bridges 6 using equation (3) produced tables 
of to and B as a function of directivity angle and jet total temperature. In addition, the results of linear 
regression provide estimates for the standard deviations cr m and a B and the correlation coefficient p m B 
between to and B. These are also tabulated. The means and the standard deviations are used to define the 
normal distribution functions when the fixed parameters are replaced with nondeterministic parameters. If 
the correlation coefficient is nonzero, then the distribution function is jointly normal. For a deterministic 
calculation, the jet operating conditions are used as input to the sJet code and the code uses interpolation 
to select the appropriate to and B values from the tables. Since the jet operating conditions do not change 
in the nondeterministic calculation, additional random inputs must be provided to the code. To simplify 
the replacement of to and B with random values, normalized random variables are provided as input and a 
linear transformation is performed. 

m = to + <r m m r m r ~ N m ( 0, 1) to ~ N m (m, a^) 

_ where => _ (8) 

B = B + a E B r B r ~N B ( 0,1) B ~ N b (B, (t 2 b ) 

The symbol N{ 0, 1) represents a normally distributed random variable with a zero mean and a variance equal 
to one. The subscripts to and B on N indicate that the random variables m r and B r are selected separately. 
The transformations in equation (8) provide the desired normal distributions for to and B where the jet 
operating conditions are used to select to, B , cr m , and a B by tabular interpolation. Up to this point, to and 
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B have been considered to be uncorrelated. In the linear regression process that produced the tables for to, 
and so forth, the correlation coefficient p m B was found to change very little and it was set at a constant value 
of 0.45 for all operating conditions and directivity angles. A constant p m B simplifies the introduction of a 
correlation between to and B into sJet. The normalized random variables m r and B r are sampled external 
to sJet with a fixed correlation and when provided to sJet the transformation maintains that correlation in 
the resulting to and B parameters. 

Following the Monte Carlo procedure, the random inputs are provided to the sJet code by DAKOTA 
using the latin hypercube sampling technique. The outputs from the sJet code are then stored and analyzed 
by DAKOTA. The first set of results compare the amplitude level given by equation (3) with fixed parameters 
determined by the linear regression to the amplitudes computed with nondeterministic parameters. Figure 2 
is a summary of results for the hot jet case at a directivity angle of 50 degrees where to = 5.75 and B = 117.76. 
The figure shows the mean levels and the 95% confidence intervals as a function of the number of random 
samples. The single point on the right is an example result after 10000 samples. In addition, the effects 
of changing er m and ag from half to twice their original tabulated value without and with correlation are 
shown. The following results are deduced from examining Figure 2: 

1. The mean value of the output distribution, here 120.5 dB, is basically unaffected by the width of the 
distributions of to and B , by whether or not there is any correlation between m and B , and by the 
number of samples. This is the expected result when the output level (in dB) is simply treated as a 
number in a linear calculation. 

2. Wider distributions in the input parameters leads to a wider distribution in the output. Again, this is 
due to the linear calculation. 

3. The width of the confidence intervals decreases as the number of samples N increases, as expected, since 
it varies as l/y/N. For the given a m and <tb and no correlation, the confidence interval is ±0.39 dB 
wide at 100 samples decreasing to ±0.12 dB wide at 1000 samples. Since this is a result of a linear 
calculation, the confidence interval is symmetric about the mean value. 

4. The presence of correlation between to and B increases the width of the output confidence interval. 
On a dB scale, this was not significant as the increase was less than 0.1 dB. 

The alternative to treating the decibel output as a number in an uncertainty calculation is to recognize 
the nonlinear nature of the power level calculation and perform the uncertainty calculation using equation 
(6). Using the same Monte Carlo approach with the normalized random input variables as above, the mean 
power levels and the 95% confidence intervals are computed before converting the results to decibels. Figure 3 
is the equivalent plot to Figure 2 for the nonlinear approach. The mean parameters for the hot jet case at a 
directivity angle of 50 degrees are to = 5.60 and E = 11.505. Using these as fixed parameters results in an 
output power level of 117.7 dB. The results for the nondeterministic calculations show the following: 

1. The mean of the output level depends on the output distribution generated by the distributions of the 
random parameters to and E. Both the size of the parameter distributions and whether the parameters 
are correlated affect the output. The fact that the mean output level increases with increasing a follows 
directly from the fact that any random change in the parameters that increases the power level of the 
output will increase the mean power in a quadratic summation of all sampled power levels. 

2. Using latin hypercube sampling, more samples results in a better definition of the output distribution 
and convergence in the statistical values computed from that distribution. Hence the mean output 
levels reach an asymptote with an increasing number of samples. 

3. The 95% confidence intervals are in general not symmetric about the mean value. Given a m , <je, and 
no correlation, the confidence interval is ±0.34 dB,— 0.37 dB at 100 samples decreasing to a confidence 
interval of ±0.11 dB,— 0.12 dB at 1000 samples. 

4. Including the effects of a correlation between parameters to and E results in an increase in the mean 
output levels and width of the 95% confidence intervals. In terms of decibels, these increases are small 
for the average correlation coefficient p rn E = 0.43, up to 0.2 dB in the mean level and up to 0.08 dB in 
the confidence interval width at 1000 samples. 
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IV. C. Many Uncertain Parameters with Unknown Distributions 

The previous model had two uncertain parameters with clearly defined distributions determined from re- 
gressions that fit model equations to measured data. We now consider the Tam model, equation (7), that 
contains 12 parameters, listed in Table 1, that each have some level of uncertainty. To determine the effects 
of these uncertainties, we have to define the probability distributions in order to replace the fixed parameters 
with distributions. For the Tam model, the detailed information required to make such distributional specifi- 
cations is not given in the references describing the development of that model. A parameter for which there 
is information available to estimate a probability distribution is the v c i parameter that is used to determine 
the convective velocity ratio U c /Uj. The data given in Harper-Bourne & Fisher 18 and Harper-Bourne, 19 
have been used to determine a fixed value for v c \ of about 0.7. That same data is used here to determine 
a mean value of 0.708 and a variance of 0.00937 for both a normal distribution and an equivalent uniform 
distribution. The uniform distribution covers a range from 0.54 to 0.876 or about ±24% from the mean 
value. Since the data were insufficient to know the exact behavior of the v c i distribution, both the normal 
and the uniform distributions were applied to replace the fixed parameter. The DAKOTA code, using the 
latin hypercube sampling technique, performed 100 samples from the v c \ parameter distribution and com- 
puted the mean and the 95% confidence interval (Cl) using the pressure-squared output distribution from 
the Tam model. The remaining parameters were fixed at the values given in Table 1. The results are shown 
in Figure 4 using the cold jet conditions at a 50 degree inlet directivity angle. For all but the lowest frequency 
results, the spectra generated using the nondeterministic v c \ parameter with either the normal distribution 
or the uniform distribution are essentially the same. Only 24 frequencies were used in the nondeterministic 
calculations, hence the mean and 95% Cl spectra are not as smooth as the fixed parameter spectrum which 
was evaluated in narrow band with a 40 Hz bandwidth. 

There are 11 other parameters in the Tam model that were identified as being replaceable by distributions. 
Without any specific information about these parameters, we estimate a distribution for each parameter 
based on the v c \ distribution results in Figure 4. First, the difference between using a normal or a uniform 
distribution was small. Second, the uniform distribution had a range based on the mean value of v c \ ± 24%. 
Observation of the v c i data indicates that this range was perhaps too wide. We reduce the range by about 
half or ±10%. Considering this to be conservative, we apply this to all the parameters listed in Table 1 
and use their listed values to define uniform distributions with a range given by their individual mean value 
±10%. A result of a 100-sample calculation is shown in Figure 5a for the hot jet case at a directivity angle 
of 50 degrees where all 12 parameters are nondeterministic and varied independently. This result shows that 
the predicted spectrum using the mean parameter values (that is the original model fixed values, shown 
by the red line) is not the same as the mean spectrum of an ensemble of spectra computed using model 
parameters whose mean values of the distributions are the same as the original fixed parameters. The 
variances or distributions of the uncertain model parameters have affected the output distribution that was 
used to compute the mean spectrum. To emphasize this, Figure 5b shows an expanded view of the main 
spectral peak where both the number of random samples and the range of the uniform distribution where 
changed. The increased number of samples provides a better definition of the output distribution. This does 
not change the output mean value significantly but does reduce the size of the confidence interval by l/\/N 
as mentioned before. Reducing the variance of the nondeterministic parameter distributions by reducing 
the range from ±10% to ±2% does show the points of the mean spectrum approaching the fixed parameter 
spectrum. 

The results in Figure 5a for the output variance from the Tam model, as expressed in the displayed 
confidence interval, when 12 parameters were varied using uniform distributions are similar to the 1 parameter 
distribution results shown in Figure 4. The confidence intervals shown in both figures are about 0.5 dB wide 
at the peak. The v c i parameter with the wider distribution has given as much output variance as all 12 
parameters give with smaller distributions. By itself, the v c i parameter with a uniform ±10% distribution 
results in a peak frequency confidence interval width of 0.12 dB. The output confidence interval for each of 
the other 11 parameters, computed separately, range from 0 to 0.68 dB in width at the peak. This indicates 
that the uncertainty in each parameter does not have the same influence on the output uncertainty and that 
computing output uncertainty due to individual parameter uncertainty neglects any interaction between 
parameters. To determine which parameters do not have significant influence and to include the effects of 
interactions between parameters, the variance-based decomposition technique is used to compute the total 
effect sensitivity measure given by equation (1). All 12 parameters have uniform distributions with a range 
of ±10% from the mean value. A result of this analysis is shown in Figure 6a for the cold jet case at a 
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directivity angle of 50 degrees. The figure shows the Sri values for each parameter added on top of each 
other as a function of the jet Strouhal number. The width of the Sti band is related to the influence of 
the parameter on the output variance. Thus, non-influential parameters have small total effect sensitivity 
measures. 

Figure 6b shows the nondeterministic parameter spectra compared to the fixed parameter spectrum and 
the measured data that corresponds to the results in Figure 6a. Comparing these two figures, we see in 
the region of the major peak that 9 of the 12 parameters have small total effect measures. The important 
parameters are those associated with the converted velocity (i> c i), the lowest mode wave number of the shock 
cell structure model (u>n), and the half-width of the similarity source (h s i). From the results in the region 
of the second peak, we include the parameter associated with the second mode wave number of the shock 
cell structure model (u> 2 i) as being an additional influential parameter. Fixing the other 8 parameters, the 
100-sample calculation with the Tam model was repeated. The resulting spectra are also shown in Figure 6b. 
The sizes of the output spectral variance are about the same between the 12-parameter and the 4-parameter 
distribution models. Thus, the sensitivity analysis has identified 8 relatively non-influential parameters in 
the prediction model. 

The same set of plots is shown in Figure 7 for the 90 degree inlet directivity angle. The peak frequency of 
the spectra have moved to a higher Strouhal number of 0.6, and the total effect sensitivity measures reflect 
this shift. The same 4 parameters as found at 50 degrees have the most influence and only the 4 parameter 
nondeterministic spectra are shown compared to the fixed parameter spectrum and the measured data. The 
50 degree and 90 degree sensitivity results are shown in Figures 8 and 9, respectively, for the hot jet case. 
These results do not show any significant changes in the influential parameters especially in the peak regions 
of the spectrum. In all these results, only the h s \ parameter has influence on amplitude uncertainty in the 
peak frequency region. The other amplitude parameters such as Cbi and abi only have their most influence on 
the output uncertainty at the higher frequencies. The parameters with consistent little influence are x c \ and 
x C 2 that are related to the jet potential core length model and v C 2 which attempts to include temperature 
effects in the model for convective velocity. 

Finally, we consider whether using nondeterministic parameters leads to better predictions. An exami- 
nation of the spectral results in Figures 6 to 9 shows that the mean spectrum computed after sampling 4 
nondeterministic parameters with uniform distributions have lowered peak levels and increased valley levels. 
At 50 degrees, the primary and secondary peaks have broadened due to the interaction of the changing wave 
numbers and the source width. This broadening reduces peak-to-valley levels. In addition, at 90 degrees the 
peak frequency of the mean spectrum shifts to a frequency lower than the fixed spectrum peak frequency and 
has better agreement with the measured peak frequency. These changes appear relatively small for the added 
complication of a nondeterministic calculation compared to the fixed parameter calculation. However, the 
total uncertainty of the Tam model output cannot be determined without varying its internal parameters. 
Consider the 4 parameters with the most influence on output spectrum uncertainty. Only the source width 
related parameter h s \ was determined from a fit to acoustic data. Its influence is both through the mode 
summation in equation (7) that determines spectral shape and directly as an amplitude parameter. The 
other parameters v c \, w\ i, and W 21 were determined in stability and waveguide models related to the jet 
flow. The true level of uncertainty in these 4 parameters is unknown and small changes in these parameters 
affect the output spectrum. If we sample these parameters using the ±10% uniform distribution and plot 
the possible spectral outputs, then after 100 samples results such as those shown in Figure 10 are obtained. 
A widely scattered set of spectra are shown where the fixed parameter spectrum is but one of the realized 
spectra. The best alternative to choosing which parameter set is best is to account for the distribution of 
the model parameters, compute the distribution of the output, and take the mean of the output distribution 
with its known level of uncertainty as your prediction. This result has not necessarily given a better pre- 
diction method, compared to any particular set of measured data, but it does provide a quantified level of 
uncertainty for the prediction and the most influential parameters were identified that may aid in developing 
an improved model. 


V. Concluding Remarks 

Engineering applications in aircraft noise prediction rely on modeling of complex physical phenomena in 
order to obtain solutions rapidly for design, development, and parametric studies. Parameters within these 
models are usually fixed in value even though there are uncertainties related to those values. It may be more 
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appropriate to allow a parameter to cover a range of values by replacing a fixed value with a distribution of 
values. A nondeterministic calculation must then be performed. In this paper, a Monte Carlo method was 
used. A random sample from one or more parameter distributions was chosen, processed by the prediction 
model, and the output value recorded. This process was repeated until a distribution of output values was 
obtained. From this output distribution, mean and confidence interval values were computed representing 
the predicted solution and its uncertainty. 

The above nondeterministic calculation process was applied to two models for the prediction of supersonic 
jet broadband shock-associated noise. The first model, called the sJet code, was developed from a study 
of far-held spectral data. It contained two nondeterministic parameters with statistical properties derived 
from performing both linear and nonlinear regression. These parameters only affect the overall level of the 
prediction and not its spectral shape. Using regression, the statistical properties of the parameters were 
obtained and used to define the parameter’s normal distribution. Random samples from these distributions 
were used to compute the output acoustic power levels. When averaging a large number of acoustic power 
levels, random higher values will raise the average level. Only the results using nondeterministic parameters 
derived from the nonlinear regression approach exhibited this behavior. 

The second model, the Tam model, was developed from a derivation of the equations of motion followed 
by using physical models to simplify the calculations. This model was used to demonstrate that the small 
uncertainty in high quality measurements of the jet operating conditions used as code input yielded small 
uncertainties in the code output. This is insufficient for demonstrating the model’s output uncertainty and 
the uncertainty of internal model parameters must be taken into account. This model contained 12 internal 
nondeterministic parameters for which estimates had to be made of their distributions in order to complete 
the analysis. External reference data were used to estimate the v c i parameter distribution. It was found that 
the output spectrum confidence interval was not changed much between using a normal distribution and a 
uniform distribution. Based on that result, all parameters were assigned a uniform distribution but with a 
smaller variance. A global sensitivity analysis performed with these 12 nondeterministic parameters identified 
8 non-influential parameters. Using the remaining 4 influential nondeterministic parameters, calculations 
were performed to judge the benefit of this approach for predictive capability. The results showed that 
the fixed parameter spectral result is only one possible solution and that the mean spectral result with its 
accompanying confidence interval computed from a distribution of output spectra would be the preferred 
result. This is the best way to quantify the uncertainty in the prediction from this model. 

Appendix A — sJet Model Spectrum Calculation 

The sJet code computes the far field radiated noise spectra from subsonic and supersonic jets. Both the 
mixing noise and the broadband shock-associated noise spectra, as appropriate, are computed as a function 
of the inlet directivity angle. For supersonic jets, the predicted radiated noise is modified depending on 
whether the jet issued from a convergent or a convergent-divergent nozzle. In addition, if screech is present, 
there is the ability to include noise amplification effects in the predicted spectra. The details of the modeling 
are found in Khavaran & Bridges. 4,6 Here, only the equations for computing the radiated broadband shock- 
associated noise spectra are given. The formulation is a modification of the model given by Deneuville 20 and 
Gliebe et al. 21 

The following parameters are defined: 

/ 3 - Shock parameter 

P = y/\M? -Mj| (9) 

where Mj is the fully expanded jet Mach number and M,j is the nozzle design Mach number with 
Mfi = 1 for a convergent nozzle. 

(3 m - Modified shock parameter for convergent-divergent nozzles, /3 m = (3 for convergent nozzles. 

For M r j > 1.01, p m depends on whether the nozzle pressure ratio Pt/Poo is greater or less than 
the critical pressure ratio. 

j (pM^ 2 \ d , -p- > pp '* , under-expanded 

Pm = < p p (10) 

I P M <i 1 , -A < p lt , over-expanded 

' MX) MX) 


NASA/TM— 2010-216801 


12 



L avg - The average spacing length between shock cells in the jet shock cell structure. 

Lavg = - 1 

where D e is the nozzle exit diameter. 

t/ c , M c - Convective velocity and Mach number 

U c = 0.7 (Uj - Uoo) 

M c = — 

Coo 

C uo - Sound level correction for convergent-divergent nozzles, M d > 1 


C un — 


-10{1.04 + tanh[4 (/3 m - 1.49)]} 


Mj > M d 


-15 jl.O + tanh 7.9 - 0.379 J j, M J < M d 


9 C - Critical angle 


7 r 

7T — tan" 


,M C < 1 


1 (yM^l), M c > 1 


For each inlet directivity angle 9e, £ = 1, 2, . . . , N ang , compute the shock-associated noise spectrum if 
< 9 C . First set the maximum sound pressure level 


SPLMAX^ = 10 nit log (3 m + Be + 42.37 + 10 log 


N s 7T / D e 

8 4V R 


i U °° , 

1 cos l 

Coo 


where N s is the number of shock cells (typically N s = 8), R is the radial distance from the nozzle to the 
microphone, and me and Be are determined by the linear regression given by equation (3). If the nozzle is 
convergent-divergent, then modify SPLMAX^ by 

SPLMAXf = SPLMAX^ + C uo Mj > M d . (17) 

The fundamental frequency of the peak of the shock-associated noise spectrum is given by 

fp,i = c —/ — jj V ' ( 18 ) 

Lavg (1 + M c cos 9 t ) ( 1 — cos 9e I 

\ Cqq / 

Finally, compute the spectrum for each frequency for k = 1, 2, . . . , K. 

t fk ( i n\ 


u * +/ 

10 log {ikt) 7 , < l 

Defp,l \ 

10 log (^) -4 , £ k e > 1 

U J 

r J J_ < 

f 101og(^/2) 2 ,£ke< 2 

’ D e f p , t 

1 10 log &</2 )- 2 ,( m >2 


SPL W = 10 log |l0 Sl/1 ° + 10 S2/10 } 
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Appendix B — Weighting in Regression Analysis 


The necessity to use weighting in the nonlinear regression formulation given in equation (6) is clue to the 
wide range of values found for yi = ( p 2 ) s h ■ The acoustic power levels considered here cover many orders of 
magnitude. 

The residual error is related to the variation in the data value yi minus the estimated value % from the 
nonlinear regression. At a particular data point after the fit, yi is a constant and the variation in the residual 
error Ci = yi — yi will follow the variation in yi. Using the analysis in Bendat & Piersol, 14 we construct an 
error band for the data value yt based on the confidence interval for an averaged spectrum. Let the following 
probability interval for the spectral density G pp 


Pr 


2 G pp <t G pp <c 


9 Gpp 

*n;l-a/2 


1 — a 


(23) 


be the definition of a (1 — a) 100% confidence interval for an averaged spectral density G pp (f) where G pp 
has a chi-squared distribution with n-degrees of freedom and Xn-a is the 100a percentage point of the % 2 
distribution. The total acoustic power is then computed as the sum of all components in the spectrum 


K 


(p 2 ) = J2 d M 

k=l 


(24) 


which is distributed as Xkxu- The confidence interval for (p 2 ) is also determined by equation (23) with 
(K x n)-degrees of freedom. The size of the variation in the data value is given by 


A (p 2 ) 


K x n 

^-Kxn;l—a /2 


( V 2 ) 


K x n 

^ K xn;a / 2 


( v 2 ). 


(25) 


Thus, if the variance of the residual error 


Var [ei] (x A(p 2 )i 


(26) 


then by using weighting and applying equation (25) 

» Var 


VVi 


Vi ~ Vi 

VVi 


Vm J 


Var[e 4 ] A (p 2 )i 

a — = constant . 

Vi [P )i 


(27) 


By using the weighting given by equation (27), the variances of the residual errors are now equal as desired 
for performing the nonlinear regression using equation (6). Furthermore, the distribution of the residual 
error is nominally x ' 2 but it approaches a normal distribution for a large number of degrees of freedom. For 
the overall power level computed from a typical power spectral density with 6000 frequency bins averaged 
400 times, the number of degrees of freedom is 2 x 400 x 6000 = 4.8 x 10 6 . 


Appendix C — Tam Model Detailed Equations 

The parametric form of the detailed model equations for the Tam broadband shock-associated noise 
prediction model, given in equation (7), are presented in this appendix. These equations follow from those 
found in Tam. 5 

c - An amplitude factor where Cn was determined by fitting the far field spectrum to measured data with 
St ref = 0.35. 

^(cfci) — £f>i x St re f (^^) 


L - The source model includes a representation of growing and decaying instability waves. L is the half-width 
of a Gaussian function that replaces the amplitude envelope of the instability wave. The coefficient h s \ was 
determined by fitting the far field spectrum to measured data. Mj^ re f = 1.67 is a reference Mach number. 


T(/r s i , x c i , x c 2^) — hsi 


X c {x cl ,x c2 ,Mj) 
x c (x cl , X c2 , M j^ref ) 


(29) 
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x c - A model for the jet potential core length. The coefficients x c i and x C 2 were determined by fitting a line 
to measured mean flow data. 


x c (x c i , X c 2 , Mj ) 


T 

Xd + X c 2 ( Mj + 1 — —A 

‘ oo 


(x c i + x C 2 Mf) Dj 


T 

Dj, :l < 1 

J ’ rp — 

OO 

Zl>1 

T 

-X oo 


(30) 


Dj - The diameter of the fully expanded jet derived by using mass conservation. 


Dj = I), 


1 + 


7 "VI fl + 1 Zllfi 


T+l 

I ^ Ty 'Ma' 1/2 

Id 


(31) 


A - An amplitude quantity that relates to the shock cell strength. The coefficients at , i and were 
determined to provide a best fit to acoustic data. 


A ( 061,062) = \ 


A 2 ( — 

1 W Mj- 

1 + 061 A“ 62 ’ M d 


> 1 


(32) 


Af 


where 


A i = 


l 1 + 061 ’ Md 

- M d\ 

l+^M* 


Mj 

> IF <! 


(33) 


/„ - The peak frequency of the spectrum associated with the u-th waveguide mode of the shock cell structure. 

, / x U c (v cl ,v c2 )k„(w 11 ,w 1 2,w2i,w„ 1 ) i 

JAVc 1 ,Vc 2 ,W 11 ,W 1 2 ,W 21 ,W v1 ) = (1 + M c (v cl , v c2 ) cos 0) (34) 

Z7T 

k„ - The wave number for the z^-th waveguide mode of the shock cell structure. The coefficients run, w 12 , 
w,' 2 i, and w„i were chosen to provide better agreement with multiple-scales shock cell calculations. 


1-1 (ton - wi 2 Mj)ki iVS , v = 1 
k„(’W 11 ,W 1 2 ,W 21 ,W„ 1 ) = { W2lk 2 ,vs ,v = 2 

rry -] kuvs , v > 2 


(35) 


where is the wave number for the vortex sheet shock cell structure and jo t „ is the zz-th zero of the zeroth 
order Bessel function. 

K.„ = (36) 

d i\M - 1 


U c - The convective velocity of large turbulent structures in the jet. The coefficients v c i and v C 2 were 
determined from stability calculations 


VclUj 


U c (v cl ,v c2 ) = 


va -v c2 \— 1 

‘ 00 


T, 


T t_ 

’ 

Tt 

Tr 


< 1 


Uj, > 1 


and the convective Mach number follows. 


M c (v c i, v c2 ) = 


Uc(v cl, V c2 ) 


(37) 


(38) 
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(a) Power spectral density at 50 deg., Mj = 1.493, 
Pt/P™ = 3.636, T t /T 0 o = 1.01 



(b) Power spectral density at 90 deg., Mj = 1.493, 
Pt/P™ = 3.636, T t /T a o = 1.01 




(c) Power spectral density at 50 deg., Mj = 1.494, (d) Power spectral density at 90 deg., Mj = 1.494, 

Pt/Poo = 3.629, Tt/Too = 2.16 Pt/Poo = 3.629, T t /T = 2.16 

Figure 1: Power spectral densities from 3 shock-associated noise prediction models compared to the shock- 
associated noise component of measured data, distance 100D e - 
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Figure 2: Mean (thick line) and 95% confidence intervals (Cl) (thin line) for sampling equation (3) with 
normally distributed m and B parameters, a range applies to both parameters where a m = 0.756 and 
<jb = 1-930. Solid lines have uncorrelated parameters. Dashed lines have p m B = 0.45. Gray area is region 
of 95% Cl when p m B = 0. Symbol on right applies to 10000 samples. 


a x 0.5, p mE = 0.00 



Figure 3: Mean (thick line) and 95% confidence intervals (Cl) (thin line) for sampling equation (6) with 
normally distributed m and E parameters, a range applies to both parameters where a m = 0.538 and 
<je = 0.177. Solid lines have uncorrelated parameters. Dashed lines have p m E = 0.43. Gray area is region 
of 95% Cl when p mE = 0. Symbol on right applies to 10000 samples. 
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Figure 4: Power spectral density comparisons at 50 degrees, distance 100D e , Mj = 1-493, Pt/Poo = 3.636, 
Tt/Too = 1.01. Effect of replacing the v c \ parameter in the Tam model with a uniform distribution ±24% 
and a normal distribution both with mean 0.708 and variance 0.00937. N = 100. Dashed lines for 95% Cl. 




(a) Uniform ±10%, N = 100 (b) Expanded view, vary N and distribution size 

Figure 5: Power spectral density comparisons showing spectra from the Tam model with 12 fixed parameters 
and spectra from the Tam model with 12 parameters replaced by uniform distributions with ±10% range 
with variance 0.00167 or with ±2% range with variance 6.68 x 10 -5 . Results at 50 degree inlet angle, distance 
100D e , Mj = 1.494, Pt/Poo = 3.629, T t /T x = 2.16. 
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Strouhal Number 


(a) Total effect sensitivity measures derived from variance-based decomposition for 
the parameters listed in Table 1. Non-influential variables have smaller widths. 
N = 10000. 



(b) Power spectral density at 50 degrees. Dashed lines for 95% CI. N = 100. 

Figure 6: Effects of reducing the number of parameters with distributions from 12 to 4 in the Tam model 
after determining variables to set constant using a variance-based decomposition sensitivity analysis. All 
parameters set with uniform distributions ±10%. Results at 50 degrees, distance 100Z? e , Mj = 1.493, 
Pt/Poc = 3.636, Tt/Too = 1.01. 
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Strouhal Number 


(a) Total effect sensitivity measures derived from variance-based decomposition for 
the parameters listed in Table 1. Non-influential variables have smaller widths. 
N = 10000. 



Figure 7: Variance-based decomposition sensitivity analysis used to reduce the number of parameters with 
distributions from 12 to 4 in the Tam model. All parameters set with uniform distributions ±10%. Results 
at 90 degrees, distance 100D e , Mj = 1.493, Pt/Poo = 3.636, Pt/Poo = 1.01. 
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Strouhal Number 


(a) Total effect sensitivity measures derived from variance-based decomposition for 
the parameters listed in Table 1. Non-influential variables have smaller widths. 
N = 10000. 



Figure 8: Variance-based decomposition sensitivity analysis used to reduce the number of parameters with 
distributions from 12 to 4 in the Tam model. All parameters set with uniform distributions ±10%. Results 
at 50 degrees, distance 100D e , Mj = 1.494, Pt/Poo = 3.629, Ft/Poo = 2.16. 
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Strouhal Number 


(a) Total effect sensitivity measures derived from variance-based decomposition for 
the parameters listed in Table 1. Non-influential variables have smaller widths. 
N = 10000. 



(b) Power spectral density at 90 degrees. Dashed lines for 95% CI. N = 100. 

Figure 9: Variance-based decomposition sensitivity analysis used to reduce the number of parameters with 
distributions from 12 to 4 in the Tam model. All parameters set with uniform distributions ±10%. Results 
at 90 degrees, distance 100D e , Mj = 1.494, Pt/Poo = 3.629, Ft/Poo = 2.16. 
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Strouhal Number 


Figure 10: Power spectral density comparisons at 50 degrees, distance 100D e , Mj = 1.494, Pt/Poo = 
3.629, Tt/Toa = 2.16. Scatter of output spectra from the Tam model computed using random samples of 4 
parameters with uniform distributions ±10% with mean output spectrum shown. N = 100. 
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